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INTRODUCTION  AND  SUMMARY 


A  calculational  parameter  study  is  presented  here 

of  the  effects  of  material  properties  upon  the  magnitude  of 

compressive  explosion-induced  residual  stresses  around  the 

cavity  produced  by  a  nuclear  explosion.  Recent  hydrofracture. 

measurements  by  W.  Ellis of  the  U.S.  Geological  Survey 

around  the  Hybla  Fair  cavity  show  that  the  rock  has  been 

noticeably  strengthened  by  the  event.  It  is  believed  that 

these  compressive  residual  stresses  are  of  primary  importance 

in  the  containment  of  the  hot  cavity  gases.  Static  two- 

dimensional  finite  element  calculations  discussed  by  Duff, 
r  ?  1 

et  al .  indicate  that  the  compressive  residual  stress 
field  greatly  increases  the  resistance  of  the  rock  to  tensile 
hydrofracture . 

This  parameter  study  involves  an  analysis  of  the  exist¬ 
ing  near  field  ground  motion  calculations  at  S3  and,  in  parti¬ 
cular,  the  late-time  calculations  made  in  connection  with  the 
teleseismic  program  of  Dr.  Cherry.  Since  only  three  calcu¬ 
lations  were  made  directly  for  this  study,  it  is  incomplete 
in  some  areas.  Some  of  the  calculations  presented  have  also 
been  discussed  in  a  previous  report  by  Rimer.  ^  However, 
attention  is  given  here  to  the  transverse  residual  stresses 
and  to  the  final  cavity  radius  rather  than  to  the  peak  radial 
residual  stresses.  Also,  the  results  of  many  more  calculations, 
not  available  previously,  are  presented  to  fill  in  some  of  the 
gaps  left  from  the  first  report. 

In  the  first  report,  an  analogy  was  made  between  the 
residual  stresses  formed  by  the  torsion  of  a  solid  steel  shaft 
or  the  bending  of  a  beam  and  the  residual  stresses  around  a 
nuclear  cavity.  It  was  hypothesized  that  residual  stresses 
are  caused  by  nonuniform  loading  into  the  plastic  regime. 

When  the  material  unloads,  stresses  result  between  elastic 
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elements,  which  cannot  recover  all  of  their  deformation. 

Any  change  in  material  property  which  tends  to  increase 
the  relative  deformation  between  elements  therefore  increases 
these  residual  stresses. 

The  results  of  this  parameter  study  tend  to  corroborate 
the  above  argument.  Some  of  these  results  are: 

(1)  Increasing  the  shear  strength  of  the  material 
increases  the  magnitude  of  the  residual  stresses. 
The  stress  magnitudes  are  particularly  sensitive 
to  the  shear  strength  at  low  pressure  (unconfined 
strength,  etc.). 

(2)  Since  increased  depth  of  burial  increases  the 
yield  strength  (for  the  modeling  used) ,  it 
also  increases  the  residual  stresses. 

(3)  Increased  shear  modulus  increases  residual 
stresses . 

(4)  A  material  which  is  stronger  in  tension  will 
have  lower  residual  stresses.  This  surprising 
result  follows  directly  from  the  above  argu¬ 
ments,  since  tensile  fracturing  increases  the 
nonuniformity  of  the  stress  distribution. 

(5)  Air-filled  voids  act  as  a  smoothing  (dis¬ 
persive)  mechanism  which  reduces  the  non¬ 
uniformity  of  the  loading  and  therefore 
reduces  residual  stresses. 

(6)  Cavity  radius  varies  inversely  as  peak  residual 
stresses  for  the  material  properties  studied. 

(7)  The  numerical  treatment  of  the  cavity  source 
(as  distinct  from  equation  of  state  of  the 
cavity)  can  influence  both  residual  stress 
and  cavity  radius.  The  cavity  treatment  used 
in  some  of  the  older  calculations  presented  in 
this  report  (and  a  recent  Pacifica  Technology 
parameter  study)  vs  more  recent  cavity  treat¬ 
ment  is  discussed  here. 
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2.  THE  EFFECT  OF  MATERIAL  PROPERTIES  ON  RESIDUAL  STRESSES 


Figure  1  shows  the  radial  and  tangential  residual 

stresses  vs  position  (compressive  stresses  are  shown  positive) 

for  a  SKIPPER  calculation  using  Husky  Ace  material  properties 

r  4  ]  . 

(see  Rimer,  et^  ajL.  )  for  a  nominal  yield  of  20  tons.  These 
profiles  correspond  to  calculation  HA14  in  Table  1  and  are 
fairly  typical  of  the  distribution  of  residual  stress  given 
by  calculations.  All  stresses  shown  in  Table  1  are  relative 
to  a  hydrostatic  overburden  pressure.  The  radial  stresses  in 
Figure  1  resist  a  cavity  pressure  (greater  than  overburden) , 
increase  with  position  to  some  peak  value  and  gradually  de¬ 
crease  to  zero.  The  peak  transverse  stress  is  generally 
greater  in  magnitude  and  occurs  closer  to  the  cavity  than 
the  radial  peak.  However,  these  transverse  stresses  decay 
rapidly  with  position,  eventually  becoming  tensile.  Thus, 
the  region  of  compressive  transverse  stresses,  which  are 
believed  to  prevent  radial  tensile  cracks  from  forming  around 
the  cavity,  is  relatively  small. 

In  this  section,  the  effect  of  calculational  variations 
in  material  properties  on  the  peak  value  of  transverse  resid¬ 
ual  stress  is  examined.  Table  1  summarizes  the  material 
properties  of  interest  and  the  resulting  peak  radial  and 
transverse  residual  stresses  as  well  as  the  cavity  radii  for 
a  nominal  20  ton  yield.  Calculations  are  identified  mnemoni- 
cally  as  follows:  Hl-Hll  and  H25-H31  refer  to  a  study  by 
J.  T.  Cherry ^  of  the  influence  of  the  material  properties 
of  a  tuff  on  teleseismic  coupling,  HR32-HR37  to  a  more  recent 
study  (see  Cherry  ^)  of  the  effect  of  water  content,  those 
labeled  HZ  to  the  same  study,  the  letter  Z  referring  to  a 
different  model  for  the  cavity  source.  Calculations  13-18 
refer  to  specific  NTS  events  detonated  in  Area  12  tuff,  HA 
to  Husky  Ace,  DS  to  Diamond  Sculls,  DQ  to  Dido  Queen,  and 
MN  to  Misty  North.  These  calculations  are  discussed  in 
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Fiqure  1  -  Radial  and  tangential  residual  stress  vs  final  position  for 
Husky  Ace  material  properties. 
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Rimer,  et  al.  YF19  refers  to  a  Yucca  Flat  dry  tuff 

simulation  discussed  in  Bache,  et  al.^  Calculations 

[3] 

S20-S24,  in  shale,  are  discussed  in  Rimer. 


2.1  CAVITY  SOURCE  MODELING 

Calculations  numbered  1  through  18,  25  through  37, 

53  through  55,  and  a  similar  parameter  study  by  Pacifica  Tech- 

r  8 1 

nology  (see  R.  Alien1  J)  used  a  cavity  source  model  which,  for 
purposes  of  this  discussion,  will  be  referred  to  as  the  ramping 
source.  This  source  introduces  a  cavity  pressure  boundary  con¬ 
dition  to  drive  the  SKIPPER  calculation.  The  cavity  pressur  is 
computed  using  an  ideal  gas  equation  of  state  given  by 


P 
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4.186  x  1019  W(Y-l) 

|-7T  R3 

3  o 


R 


3y 


o 


R  (1) 


(dynes/cm  ) 


(1) 


where 

W  =  device  yield  (Kilotons) 

R  =  initial  cavity  radius  (cm) 

o 

R (1 )  =  cavity  radius  (cm)  at  time  t 
y  =  ratio  of  specific  heats. 


The  cavity  pressure  is  "ramped,"  i.e.,  restrained  to 
build  up  slowly  for  the  first  50  to  100  cycles  by  modifying 
the  above  equation  as  follows 

P  =  P  P  <  P  (2) 

c  t  c 

o 

where  t  is  the  time  it  would  take  for  a  wave  propagating  at 
o 

the  sound  speed  of  the  cavity  material  at  time  t  to  propagate 
through  4  zones.  Ramping  is  discontinued  and  Equation  (1)  is 


used  when  t  remains  greater  than  t  for  20  consecutive 
cycles.  This  is,  of  course,  artificial. 

The  original  purpose  of  this  ramping  source  model  is 
related  to  the  Tabular  Arrays  of  Mixtures  Equation  of  State 
(TAMEOS) ,  described  later  in  this  report.  The  original 
TAMEOS  table  used  in  the  calculations  only  reached  to  a  maxi¬ 
mum  pressure  of  600  kbars.  Ramping  was  used  to  avoid  having 
the  pressure  in  the  first  grid  zone  outside  the  cavity  exceed 
the  bounds  of  the  table.  When  tables  were  developed  reaching 
to  higher  maximum  pressures,  ramping  was  no  longer  needed. 
However,  in  order  to  compare  older  with  newer  calculations, 
the  ramping  model  was  retained  far  longer  than  necessary. 

The  ramping  source  model  has  many  disadvantages,  the 

most  important  being  that  it  is  nonconservative  during  the 

ramping  phase.  The  amount  of  ramping  and  therefore  the  energy 

loss  is  dependent  upon  the  zoning  used  in  the  calculation 

(tQ  depends  on  zoning) .  Energy  losses  of  10  to  15  percent 

have  been  observed  in  the  calculations  presented  in  this 

report.  However,  losses  of  40  to  50  percent  have  been  seen 

in  calculations  in  sand  (40  percent  air-filled  voids), 
f  9 1 

R.  B]ork  considers  the  energy  remaining  in  the  calculation 

after  the  ramping  phase  to  be  the  actual  device  yield  of  the 

calculation.  This  is  probably  the  best  approach  to  use,  but 

does  not  eliminate  the  difficulty  associated  with  ramping. 

Comparisons  presented  in  this  section  indicate  that  cavity 

radii  obtained  with  and  without  ramping  differ  by  far  more 

1/3 

than  the  energy  lost  (using  W  scaling) .  This  is  also  true 
of  reduced  displacement  potentials  (RDP) ,  measures  of  the 
coupling  efficiency  of  elastic  waves  from  a  device  into  the 
ground.  The  RDP  scales  with  yield  and  increases  by  approxi¬ 
mately  35  percent  when  ramping  is  not  used. 


As  a  first  step  toward  improving  the  cavity  represen¬ 
tation,  the  so-called  "bubble"  model  (see  Schroeder ^ )  was 


incorporated  into  SKIPPER.  This  model  assumes  that  70  metric 
tons  of  rock  at  ambient  density  is  vaporized  for  each  Kiloton 
of  nuclear  explosive  in  order  to  compute  an  initial  cavity 
radius.  The  device  energy  is  loaded  into  the  cavity  which  is 
modeled  using  one  grid  zone  with  an  ideal  gas  equation  of 
state  with  constant  y  of  1.5.  This  procedure  results  in 
energy  conservation  to  better  than  one  percent  in  all  calcu¬ 
lations.  In  this  report,  calculations  19-24  and  38-51  were 
made  using  this  cavity  model.  More  recent  work  on  a  cavity 
equation  of  state  by  D.  Laird  indicates  that  y  varies 

from  approximately  1.55  at  zero  time  to  1.1  at  later  times. 
Since  the  influence  of  the  cavity  at  early  times  is  far  more 
important  to  the  wave  propagation,  the  choice  of  y  =  1.5 
seems  reasonable  as  a  first  approximation. 

Calculations  1  through  12  and  25  through  37  and  HA52 
for  the  ramping  source  model  used  y  =  4/3  for  the  source 
material,  while  13  through  18  used  y  =  1.5.  All  used  an 
arbitrary  initial  cavity  radius  which  is  equivalent  to  vapor¬ 
izing  approximately  35  tons  per  Kiloton  of  explosive.  Com¬ 
parisons  between  the  bubble  model  and  the  ramping  source  may 
be  made  from  Table  1.  Three  comparisons  are  reproduced  in 
Table  2,  each  pair  for  identical  material  properties  outside 
the  cavity  source.  In  all  cases,  the  bubble  model  used 
y  =  1.5  in  the  cavity  (y  for  the  ramping  model  varied).  The 
results  of  calculations  HA13  for  y  =  1.5  and  HA52  for 
y  =  4/3,  which  have  identical  material  properties,  show  that 
at  least  for  the  ramping  source  model,  the  value  of  y  has 
negligible  influence  on  residual  stresses.  The  larger  y 
did  increase  the  cavity  radius  by  about  2  percent.  Due  to 
the  apparent  small  effect  of  changing  y,  the  variation  in 
y  in  Table  2  will  be  ignored. 

Comparisons  1  and  2  indicate  only  slight  differences  in 
residual  stresses.  Cavity  radii  for  the  bubble  model  are 
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approximately  12  percent  greater  in  both  cases-  A  similar 
increase  in  cavity  radius  is  seen  in  comparison  3.  However, 
unlike  the  first  two  comparisons,  residual  stresses  show 
significant  increases  for  the  bubble  calculation. 

Summarizing,  the  bubble  cavity  source  model  is  an 
improvement  over  the  ramping  source  model  which  is  noncon¬ 
servative.  The  bubble  model  leads  to  larger  cavity  radii; 
however,  further  work  is  needed  to  determine  the  effect  of 
cavity  modeling  on  peak  residual  stresses.  For  the  parameter 
study  to  be  presented  in  this  report,  comparisons  between 
calculations  having  different  cavity  source  models  will  be 
avoided. 

2.2  WATER  CONTENT  AND  AIR-FILLED  VOIDS 

All  calculations  presented  except  S20  through  S24 

utilize  the  TAMEOS  tabular  equation  of  state  for  tuff.  A 

TAMEOS  table  is  generated  for  a  given  water  content  (mass 

fraction  of  water)  assuming  zero  air-filled  porosity.  The 

table  describes  all  relevant  states  of  the  mixture  of  dry 

(crystalline)  tuff  of  grain  density  2.4,  represented  by  an 

equation  of  state  of  the  Mie-Gruneisen  type  and  water,  repre- 

f  121 

sented  by  an  equation  of  state  developed  by  Bjork.  The 

TAMEOS  scheme,  described  in  Riney,  et_  al. ,  assumes  pressure 

equilibrium  between  the  dry  tuff  and  water. 

The  air-filled  voids  are  included  through  the  S3  con¬ 
stitutive  equation  for  porous  materials  (P-a  model)  shown 
graphically  in  Figure  2.  The  pressure  of  the  porous  material 
is  described  by 
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SKIPPER 


(sjeqq)  ajnssaJd 
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Figure  2  -  Loadina  and  release  p-V  curves  for  partiallv 
saturated  tuff  (f  =  0.]7,  4>  =  0.05)  [cherry 


where 


v  =  specific  volume  of  the  material 
e  =  specific  internal  energy 
P  =  pressure  obtained  from  the  TAMEOS  table 
a  =  distension  ration  defined  by 


where  v  is  the  specific  volume  of  the  material  with  zero 
air-filled  voids.  The  distension  ratio  is  required  to  decrease 
from  an  initial  value  at  zero  pressure  down  to  1.0  as  the 
pressure  increases  to  P  ,  the  crush  pressure  at  which  all 
air-filled  porosity  is  irreversibly  removed.  The  pressure 
limit  for  the  completely  reversible  portion  of  the  void  col¬ 
lapse  ?e,  locates  the  boundary  between  two  functions  a(v), 
elastic  and  plastic,  which  together  define  the  crush  curve. 
Unloading  from  the  plastic  portion  of  the  curve  allows  for 
partial  pore  recovery  in  order  to  provide  continuity  at  P  . 

Thus,  the  dotted  curves  in  Figure  2  are  not  parallel. 

Figure  3  shows  the  effect  of  water  content  on  the  peak 
value  of  transverse  residual  stress.  Each  curve  has  all 
material  properties  except  water  content  held  constant,  the 
numbers  on  a  curve  representing  the  corresponding  calculation 
numbers  in  Table  1.  Where  a  curve  consists  of  only  two  data 
points,  they  are  joined  by  a  straight  line.  This  is  not 
intended  to  imply  a  straight  line  relationship.  In  fact,  curve 
2,  connecting  the  results  of  calculation  H10  (for  zero  water 
content)  with  H4  for  17  percent  water  by  weight,  appears  to  give 
a  totally  erroneous  impression  of  the  variation  of  residual  stress 
with  water  content. 
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Figure 


The  data  points  for  zero  water  content  on  curves  1  and  3 
(connected  to  the  rest  of  the  curves  by  dotted  lines)  do  not 
appear  to  be  consistent  with  the  remainder  of  the  parameter  study 
which  shows  a  general  increase  in  peak  residual  stress  with  water 
content  for  low  water  contents.  Numerically,  the  calculations 
for  zero  water  content  differ  from  the  other  calculations  in 
that  the  TAMEOS  scheme  is  not  used.  Mixing  water  with  the  dry 
tuff  using  TAMEOS  generates  an  equation  of  state  having  releases 
very  different  from  the  dry  tuff  alone.  It  is  not  clear  at  this 
time  why  this  influences  the  residual  stress  variation  at  low 
water  contents. 

Curve  3  of  Figure  3  shows  very  little  change  in  peak 
residual  stress  for  mass  fractions  of  water  above  7  percent. 
However,  one  very  important  effect  of  water  content  was  not 
included  in  this  parameter  study;  the  decrease  in  the  shear 
strength  of  the  rock  due  to  increased  Water  content.  Results 
will  be  presented  later  in  this  report  that  show  that  decreas¬ 
ing  strength  decreases  the  peak  residual  stress  dramatically. 
Thus,  one  must  draw  conclusions  from  the  results  of  these  cal¬ 
culations  with  care.  The  influence  of  water  on  equation  of 
state  has  been  included.  The  influence  of  water  on  rock 
strength  has  not.  It  may  well  be  that  the  net  effect  of 
increasing  water  content  will  be  a  decrease  in  residual  stress 
field  around  the  cavity. 

Figure  4  shows  the  dependence  of  residual  stresses  on 
the  volume  fraction  of  air-filled  voids.  The  curve  summarizes 
calculations  HI,  H6  and  H8  and  indicates  that  the  peak  residual 
stresses  are  greater  for  materials  with  less  voids.  The 
crushing  up  of  these  voids  acts  as  a  smoothing  (dispersive) 
mechanism  which  may  reduce  the  nonuniformity  of  the  plastic 
loading.  An  argument  has  already  been  presented  that  it  is 
this  nonuniformity  which  causes  residual  stresses. 
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Figure  4  -  Peak  transverse  residual  stress  vss  air 
filled  voids. 


Calculations  13  to  17  differ  only  slightly  in  air-filled 
voids.  In  order  to  compare  other  features  of  these  calcula¬ 
tions.  we  will  assume  that  these  calculations  do  not  differ 
in  this  property.  Figure  4  indicates  that  the  difference 
between  1.2  and  1.6  percent  air-filled  voids  is  only  4  percent 
in  transverse  residual  stress.  Similarly,  the  difference  in 
water  contents  among  these  calculations  is  only  a  small  effect. 

The  comparison  of  calculations  Hll  and  H27  shown  in 
Figure  5  indicates  that  a  higher  crush  pressure  results  in 
a  larger  residual  stress.  It  reduces  the  smoothing  effect 
by  steepening  the  crush  curve  and  by  moving  the  crush  curve 
into  the  higher  stress  regime.  In  a  similar  way,  the  increase 
in  residual  stresses  due  to  the  increase  in  elastic  pressure 
shown  in  Figure  5  may  be  explained. 

For  all  the  nonsaturated  materials  studied,  the  maximum 
residual  stresses  were  observed  near  the  radial  position  cor¬ 
responding  to  the  end  of  the  fully  crushed  region  (where  peak 
radial  stress  is  equal  to  the  crush  pressure).  Material 
loaded  to  above  Pc  unloads  along  its  isentrope  independent 
of  the  shape  of  the  crush  curve.  However,  material  which  was 
never  fully  crushed  has  an  unloading  path  which  depends  on 
the  shape  of  the  crush  curve  and  on  the  amount  of  loading 
it  has  been  subjected  to.  This  nonuniformity  between  load- 
unload  paths  for  the  materials  near  Pc  leads  to  large  final 
state  gradients  and  is  probably  why  the  peak  residual  stress 
is  located  at  the  end  of  the  fully  crushed  region. 

2.3  SHEAR  MODULUS 

Curves  1,  2  and  3  of  Figure  6  indicate  that  increasing 
the  shear  modulus  results  in  an  increase  in  peak  transverse 
residual  stress.  Curve  4,  the  results  of  calculations  HZ41, 
HZ49  and  HZ50,  however,  show  a  very  slight  decrease  in  resi¬ 
dual  stress.  The  effect  of  shear  modulus  appears  small,  but 
is  very  difficult  to  decouple  from  the  effect  of  varying  zero 
pressure  bulk  modulus  and  sound  speed  which  are  related  by 
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transverse  residual  stres 


Elastic  pressure,  (kbars) 


Peak  transverse  residual 


o  C2  =  K  +  4/3  G 
o  o  o 


where 


p  =  density 
C  =  sound  speed 
K  =  bulk  modulus 
G  =  shear  modulus 

All  calculations  shown  use  a  constant  shear  modulus 
to  relate  deviatoric  stress  history  to  strain  rates.  Thus, 
the  major  influence  of  increased  shear  modulus  is,  for  the 
same  strain  history,  to  reach  the  von  Mises  failure  surface, 
i.e.,  to  behave  plastically,  more  quickly.  Residual  stresses 
will  be  shown  to  be  coupled  very  closely  to  the  plastic 
behavior  of  the  rock  and  to  the  failure  surface. 


2.4  FAILURE  SURFACE 

The  material  strength  model  used  in  SKIPPER  requires 
that  the  principal  stress  be  within  the  von  Mises  yield 
surface.  This  is  equivalent,  for  spherical  symmetry,  to 
the  radial  deviatoric  stress  SD  satisfying  the  following 
inequality 

I SD |  <  2/3  Y  (5) 

where  Y  is  the  yield  strength  (maximum  stress  difference) 
equivalent  to  twice  the  maximum  shear  strength.  The  non- 
associated  flow  rule  is  used  for  all  calculations.  Thus, 
whenever  SD,  as  calculated  from  the  strain  rates,  becomes 
greater  in  absolute  value  than  2/3  Y,  its  magnitude  is 
decreased  to  2/3  Y.  The  yield  strength  Y  ,  considered  to 
depend  both  on  pressure  and  energy,  is  given  by 
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Y  =  0  ,  e  .>  e  (6c) 

m 

where 

Yq  =  cohesive  strength 

Y  =  maximum  increase  in  Y  due  to  pressure 
m 

P  =  pressure  at  which  the  maximum  stress  difference 

Ll  under  triaxial  compression  (Y  +  Y  )  is  attained 

o  m 

em  =  melt  energy 

Here,  P  includes  a  scalar  overburden  pressure. 

A  comparison  between  calculations  HR37  and  HR33  shows 
that  increasing  the  melt  energy  by  over  a  factor  of  2  decreases 
the  peak  transverse  residual  stress  by  only  about  3  percent. 
Thus,  em  does  not  appear  to  be  an  important  parameter. 

Figure  7  shows  four  comparisons,  all  of  which  indicate 

that  increasing  the  maximum  yield  strength  (Yq  +  Y^)  increases 

the  magnitude  of  the  residual  stresses.  For  the  analogy  of 

torsion  of  an  ideal  plastic  cylinder,  the  magnitude  of  the 

peak  residual  stress  can  be  shown  analytically  to  be  linearly 

r  3 1 

proportional  to  the  yield  stress  (see  Rimer1  ).  Thus,  the 
strong  dependence  between  strength  and  residual  stresses 
around  a  nuclear  cavity  is  not  surprising. 

The  curves  in  Figure  7  are  for  different  failure  sur¬ 
faces.  Some  have  nonzero  cohesive  strengths,  while  others 
do  not.  Figure  8  shows  the  effect  of  varying  Y  while 
maintaining  the  maximum  yield  strength  (Yq  +  Y^)  constant  at 
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Figure  8  -  Peak  transverse  residual  stress  vs 
cohesive  strength  (Y  ) . 
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1  kbar.  Increasing  the  proportion  of  Y  in  the  yield  strength 

increases  the  residual  stresses  since  for  any  given  pressure 

less  than  P  ,  Y  will  be  increased.  The  limiting  case  of  no 
m 

pressure  dependence,  Y  =  Yq  (calculation  HZ45)  is  shown  for 
comparison . 

2.5  DEPTH  OF  BURIAL 

The  depth  of  burial  of  the  device  is  represented  in 
this  study  by  an  overburden  pressure  which  is  added  to  the 
pressure  computed  from  TAMEOS.  One  effect  of  overburden  is 
to  increase  the  work  done  term  and  therefore  to  increase  the 
internal  energy  further  during  compression.  The  most  important 
effect  of  overburden  on  residual  stress  is  its  influence  on  the 
failure  surface  given  in  Equation  (6) .  Overburden  is  added 
onto  the  pressure,  increasing  the  yield  strength.  Figure  9 
shows  the  resulting  increase  in  peak  residual  stress  with  over¬ 
burden  . 

2.6  TENSION  CRACKS 

Calculations  S20  through  S24,  for  shale,  use  the  modi¬ 
fy] 

fied  Tillotson  equation  of  state  discussed  in  Rimer.  These 

calculations  allow  the  material  to  crack  in  tension  whenever  a 
principal  stress  becomes  greater  than  some  maximum  tensile 
stress.  The  crack  is  allowed  to  grow  and  subsequently  close 
when  the  stress  distribution  becomes  compressive.  Figure  10 
indicates  that  the  higher  the  allowable  stress  in  tension  is 
made,  i.e.,  the  stronger  the  material  in  tension,  the  weaker 
the  residual  stress  field  around  the  cavity  becomes.  It 
appears  that  whenever  tension  cracks  open  in  the  rock  and 
subsequently  close,  the  residual  stress  membrane  around  the 
cavity  is  stronger  than  if  no  cracking  occurred.  Calculation 
S21  represents  the  effect  of  no  cracking  and  corresponds  to  a 
tensile  strength  greater  than  or  equal  to  27  bars. 


25 


Peak  tra: 


The  above  results  may  be  explained  in  terms  of  distri¬ 
bution  of  stress  around  the  cavity.  When  tensile  cracks  are 
opened,  the  cavity  region  where  residual  stresses  form  is 
essentially  decoupled  from  the  propagating  stress  waves.  Thus, 
cracking  increases  the  nonuniformity  of  the  stress  distribution 
in  this  region.  Since  rock  has  been  loaded  into  the  plastic 
regime,  on  unloading,  greater  residual  stresses  should  be 
present . 
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3.  THE  EFFECT  OF  MATERIAL  PROPERTIES  ON  CAVITY  RADIUS 

Fjnal  cavity  radii  are  presented  in  Table  1  for  a 
nominal  device  yield  of  20  tons.  From  these,  cavity  radii 
mav  be  obtained  for  other  yields  using  scaling.  The 

effect  of  material  Properties  on  cavity  radius  is  studied 
ne-e  *nd  a  correlation  made  with  residual  stresses.  The 
effect  o*  cavity  source  modeling  has  been  discussed  in 
Chapter  2. 

5.1  WATER  CONTENT,  AIR-FILLED  VOIDS  AND  SHEAR  MODULUS 

Figure  11  indicates  that  increasing  the  water  content 
decreases  the  final  cavity  radius.  An  effective  stress  law 
would  be  expected  to  qive  an  increase  rather  than  the  decrease 
seen  in  Figure  11.  However,  the  decrease  in  strength  due  to 
increased  water  content  is  not  included  in  the  calculations 
used  for  this  study.  Therefore,  the  results  shown  must  be 
duo  to  the  complicated  releases  from  high  pressure  states 
given  by  TAMEOS  for  water-rock  mixtures. 

Increasing  air-filled  voids  results  in  the  increase 
ir,  cavity  radius  shown  in  Figure  12.  When  the  voids  are 
crushed  up,  the  cavity  has  more  room  to  expand.  Increasing 
the  crush  pressure  pc  results  in  less  void  crushup.  There¬ 
fore,  as  shown  in  Figure  13,  cavity  radii  are  smaller.  Simi¬ 
larly,  a  higher  elastic  pressure  pe  results  in  less  permanent 
crushup  and  smaller  cavity  radius. 

Figure  14  indicates  that  higher  shear  modulus  materials 
give  smaller  cavity  radii. 

3.2  FAILURE  SURFACE,  DEPTH  OF  BURIAL  AND  TENSION  CRACKS 

The  modeling  of  the  yield  surface  was  discussed  in 
Section  2.4.  Figure  15  shows  that  yield  strength  is  the 
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Figure  11  -  Final  cavity  radius  vs  mass 
fraction  of  water. 
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Figure  14  -  Final  cavity  radius  vs  shear  modulus. 
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Cavity  radius  (meters) 


material  property  of  greatest  influence  on  cavity  radius. 

Strong  materials  give  small  cavities.  This  is  in  agreement 

[14] 

with  earlier  work  by  Allen  and  Duff.  Figure  16  indicates 

that  the  cohesive  strength  Yq  is  very  influential  in  deter¬ 
mining  final  cavity  radius.  An  argument  may  be  made  that  at 
early  times,  the  cavity  pushes  out  so  strongly  that  it  does 
not  even  notice  the  strength  of  rock  around  it.  However, 
when  the  surrounding  rock  has  unloaded  to  lower  stresses,  the 
strength  at  low  pressure  is  predominant  in  slowing  down  and 
finally  stopping  cavity  growth.  If  the  rock  has  negligible 
strength,  the  cavity  will  continue  expanding  forever.  In 
all  calculations,  it  is  noticed  that  the  cavity  expands 
beyond  its  final  radius  and  then  reverses  its  motion.  It  is 
during  this  time  period  that  the  residual  stresses  reach  their 
final  magnitudes. 

Increased  depth  of  burial,  for  the  modeling  used,  re¬ 
sults  in  higher  yield  strength  at  lower  dynamic  pressures. 
Therefore,  increasing  overburden  pressure  results  in  decreases 

in  cavity  radii.  These  are  quite  close  to  those  predicted  by 
-1/4 

the  CEP  h  relationship  as  shown  in  Figure  17.  Actually, 

the  calculated  cavity  radius  does  not  decrease  quite  as  fast 

-1/5 

as  the  CEP  formula  suggests.  An  h  dependence  is  in 

better  agreement  with  the  calculational  results. 

Figure  18  shows  no  clear  trend  between  allowable  tension 
stress  and  cavity  radius. 

3.3  CAVITY  RADIUS  AND  PEAK  RESIDUAL  STRESS 

Figure  19  is  a  plot  of  peak  transverse  residual  stress 
vs  cavity  radius;  each  point  in  the  plot  represents  one  of  the 
calculations  of  Table  1.  The  trend  is  clear.  A  larger  cavity 
radius  indicates  a  smaller  residual  stress.  When  variations 
in  an  individual  material  property  are  examined,  this  trend  is 
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Fiqure  17  -  Final  cavity  radius  vs  overburden.  (Curves 


showing  variations  with  depth  of  burial  of 
h"'*'//^  and  are  included  for  reference.) 
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Fiqure  19  -  Peak  transverse  residual  stress  vs  final 
cavity  radius  for  calculations  of  Table  1 


seen  for  all  properties  with  the  exception  of  allowable  tensile 
strength  where  no  clear  trend  is  evident. 

It  has  been  noted  that  the  residual  stresses  reach  their 
final  values  at  about  the  time  of  cavity  formation.  The  cavity 
generally  overexpands  and  rebounds  inward.  At  this  time,  the 
rock  is  squeezed  together,  due  to  convergence  effects,  giving 
rise  to  the  large  hoop  stresses  (transverse  residual  stresses) . 
This  convergence  effect  is  greatest  for  smaller  cavity  radii. 
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